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The renormalization-group (RG) approach proposed earlier by Shankar for interacting spinless fermions at T = 
is extended to the case ol non-zero temperature and spin. We study a model with S'L^AQ-invariant short-range 
effective interaction and rotationally invariant Fermi surface in two and three dimensions. We show that the 
Landau interaction function of the Fermi liquid, constructed from the bare parameters of the low-energy effective 
action, is RG invariant. On the other hand, the physical forward scattering vertex is found as a stable fixed point 
of the RG flow. We demonstrate that in d = 2 and 3, the RG approach to this model is equivalent to Landau's 
mean-field treatment of the Fermi liquid. We discuss subtleties associated with the symmetry properties of the 
scattering amplitude, the Landau function and the low-energy effective action. Applying the RG to response 
functions, we find the compressibility and the spin susceptibility as fixed points. 
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I. INTRODUCTION 

The necessity to better understand the Physics of 
strongly correlated fermions, like in high-T c supercon- 
ductors or in the fractional quantum Hall effect, inspired 
a lot of theoretical efforts to explain the observable de- 
viations from Fermi liquid (FL) behavior, and to clarify 
the foundations of the Landau Fermi liquid theory itself. 
Much work has been done to elucidate the microscopic 
basis olihe FL theory in the framework of quantum field 
theoryda. However, the development of Renormaliza- 
tion Group (RG) methods to the problem of interacting 
fermions in dimension d is relatively new. 

So far, the most successful applications of RG meth- 
ods to interacting fermions have been achieved in the 
one-dimensional case, where it was known from exact 
solutions that non FL phases exist (e.g. the Luttinger 
liquid). For a review on ID systems see Ref. [| Later, 
Bourbonnais and CaronP made an extensive RG study of 
one-dimensional and quasi-one-dimensional fermion sys- 
tems at finite temperature, in which a transition from 
Luttinger liquid to FL behavior is revealed. The de- 
velopment of a RG theory for interacting fermions in 
isotronkjjSMStems of dimensions greater than one is more 
recentu'EffilQ. Such a generic theory, when completed, 
should include the Landau FL theory as one of the fixed 
points of the RG equations. As discussed by Shankar in 
a very pedagogical paper p the RG treatment of fermions 
is a much more complicated issue than the analogous 
procedure for critical phenomena, because of two crucial 
points: Firstly, the existence of a Fermi surface: the low- 
energy modes lie in the vicinity of a continuous geomet- 
rical object (the Fermi surface) and not only around iso- 
lated points, like the origin of phase space in the bosonic 
case. This introduces additional phase space constraints 
on the modes to be integrated out, quite a problem for an 
arbitrary Fermi surface. Moreover, the Fermi surface it- 



self is a relevant parameter of theory, and its shape should 
renormalize under the RG procedure, except in the ro- 
tationally symmetric case. Secondly, instead of studying 
the flow of one or a few coupling constants like in criti- 
cal phenomena, one has to set up RG flow equations for 
coupling functions, defined on the Fermi surface. Notice 
that the purely ID interacting fermion system can be 
considered as a degenerate special case where the Fermi 
surface is reduced to a set of two points, with a finite 
number of coupling constants. It is therefore closer to 
the usual applications of the renormalization group (how- 
ever, in the quasi-one-dimensional approach of Ref. 0, one 
already witnesses the appearance of new coupling con- 
stants related to interchain hopping and one may follow 
the change in shape of the (open) Fermi surface under 
RG flow). 

In this paper we extend Shankar's approachtl, devel- 
oped for spinless fermions at zero temperature. The main 
goal of our analysis is to incorporate spin and a finite tem- 
perature. For the sake of generality, we study a system of 
A^-component fermions with an SU (A^)-invariant short- 
range interaction. We perform the analysis for a circular 
Fermi surface in dimension two, spherical in dimension 
three. 

Shankar's analysis of coupling functions for the quar- 
tic interaction shows that only two such functions survive 
under the RG flow: the function $, which couples two 
incoming and two outgoing particles with the same pairs 
of momenta lying on the Fermi surface (Landau channel) , 
and the function V, which couples incoming and outgo- 
ing particles with opposite momenta lying on the Fermi 
surface (Bardeen-Cooper-Schrieffer (BCS) channel). 

We find that the coupling function V is marginally ir- 
relevant if its bare value satisfies Landau's conditions for 
the stability of the Fermi liquid against Cooper pairing at 
arbitrary momentum; otherwise a BCS instability occurs 
at finite temperature. In the Landau channel the distinc- 
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tion between the Landau interaction function <t and the 
physical forward scattering vertex T appears naturally in 
the finite temperature RG formalism. The function $ is 
RG invariant at non-zero temperatures, while the vertex 
f is attracted towards a non trivial stable fixed point. In 
order for the Fermi liquid to be stable up to zero tempera- 
ture, we find that the components of the Landau function 
must obey conditions-which coincide with Pomeranchuk's 
stability conditions£3 

The paper is organized as follows: In Section II we de- 
scribe the low-energy effective action which is the start- 
ing point of our RG analysis. In Section III we study the 
two-dimensional case. This section contains a description 
of the key points of our analysis and most of the techni- 
cal details. There, we obtain the RG equations for the 
vertices and their fixed points. In Section IV we apply 
the same analysis to d = 3. In Section V we calculate 
the compressibility and the spin susceptibility in this RG 
framework. 



II. THE LOW-ENERGY EFFECTIVE ACTION 

We treat the problem of interacting fermions at Anile, 
temperature in the standard path integral formalism^ 
using Grassmann variables for Fermi fields. The partition 
function is given by the path integral 

Z = J V^VtP e So+s ' nt (1) 

wherein the free part of the action is 

So= f &,(l)[»Wi + M-e(Ki)]^a(l) (2) 
J(i) 

We introduced the following notation: 

l - um 

(i) = (Ki.wO (3b) 

where (3 is the inverse temperature, fi the chemical po- 
tential, oji the fermion Matsubara frequencies and ifi a (i) 
an V-component Grassmann field with a 'flavor' index a. 
Summation over repeated indices is implicit throughout 
this paper. We set fee = 1 and % = 1. The physically 
interesting case is of course N = 2, but the generalization 
to N ^= 2 is simple, and it incorporates automatically the 
simpler case of spinless fermions (N = 1). 

The general S'C/(A r )-invariant quartic interaction may 
be written as follows: 

Sint = -t / ^ Q (1)#(2)V 7 (3)^(4) 

4 ^(1,2,3,4) 

x[Tf (1, 2; 3, 4)<^ +1 )(l + 2-3-4) (4) 

Here the conservation of energy and momentum is en- 
forced by the symbolic delta function 



S (d+1 \l + 2 - 3 - 4) = f3(2ir) d x 

<5(Ki + K 2 - K 3 - K 4 )A(wi + uj 2 - to 3 - w 4 ) (5) 

where the discrete delta function A is equal to 1 if its 
argument is zero, and equal to zero otherwise. 

It can be shown in group theory (see, for instance, 
Ref. |l2|) that a representation of the potential as 

N 2 -l 

U$ = UxSatS^ + £/ 2 ]T W 7 (6) 

a=l 

supplies us with the most general SU (TV)-invariant form 
with two independent scalar functions U\ and Ui. The 
N 2 — 1 Hermitian traceless matrices A a are the generators 
of SU(N). We need not write down explicit expressions 
of the commutations or other relations for those matrices. 
The only identity used in the following is 

X! A a/3 A 7<5 = 2 ( SasSp-y ~ J (7) 

o=l ^ ' 

In the SU(2) case, the three generators A a are the usual 
Pauli matrices and the identity (^) reduces to a well- 
known relation involving these matrices. Using the rela- 
tion (0) , one readily checks that the form (^|) of the inter- 
action is indeed SU(N) invariant. The fact that Eq. (||) 
is the most general SU (N) invariant may be verified by 
counting the number of singlets in the tensor product 
N <£> N <g> N <E> N , wherein N stands for the fundamen- 
tal representation of SU (N) (acting on the iV-component 
field VO and N for its conjugate. This number is indeed 
two, meaning that only two SU (N )-invariant scalars may 
be constructed in this way. 

One should bear in mind the difference between the 
5L r (A r )-invariant interaction considered here and a rota- 
tion invariant interaction for particles of spin s: In the ab- 
sence of symmetry-breaking interactions (e.g. spin-orbit, 
dipole-dipole, or external fields) there will be rotational 
invariance in spin space, but the corresponding symme- 
try operations are still obtained from the three genera- 
tors of SU(2), although in a (2s + l)-dimensional repre- 
sentation. The 5'L r (7V)-invariance imposed here is more 
stringent. The particular form of Eqs. (|]|7|), and conse- 
quently of Eq. (|ll]), for any N, is the artifact of such a 
symmetry enlargement. A generic, 5J7(2)-invariant in- 
teraction with spin-s fermions would contain a greater 
variety of terms. Accordingly, the different components 
of ip are called 'flavors' if N = 2s + 1 > 2, in order to 
avoid misunderstandings. 

The flavor dependence of the potential (||) may be fac- 
torized and expressed via two independent flavor opera- 
tors. It is convenient to introduce two operators / and 
T, respectively symmetric and antisymmetric, as follows: 

I"s = SasSpj + Sa^Sps (8a) 
T"g = SasSp^ - SaySps (8b) 
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These operators satisfy the properties 



A < K F 



(13) 



/ 7 f = J# = J# (9a) 



-t: 



a/3 
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rpO.fi a^/3a 

J 7 <5 — ~ 1 1 8 

and the convolution relations 

rap jv[i _ ra/3 
1 iiu 1 1 8 — Zll 1 S 

J F J 7 i — ZJ 7 <5 

rpafi TVIl _ n 

1 t iv 1 1 S — U 

TCtfj, rvfi _ N + 3 af3 _ N + 1 a p 
fv fj,S 2 7* 2 ' 

TO^irpVfi _ N - I a p N + 1 gff 

fiy ^lS 2 7^~ 2 7 d 



rpOllrpVfi 

fv fj.8 



N - 1 



a/3 



- tU/_ 

~ J 7<5 



AT-3. 



a/3 



" J 7(5 



(9b) 

(10a) 
(10b) 
(10c) 

(lOd) 
(lOe) 
(lOf) 



Instead of Eq. (g|), we may decompose the potential as 
follows: 



jja{3 _ jjAja/3 ^ jjSrpa/3 



7<5 



7<5 
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(ii) 



where the functions U s and U have the symmetry prop- 
erties 

C/ A (l,2;3,4) = -£/ A (2,l;3,4) = -U A (1, 2; 4, 3) (12a) 
t/ s (l,2;3,4) = C/ S (2,l;3,4) = £/ s (l,2;4,3) (12b) 

The general form (Q) of the interaction allows us to eas- 
ily recover various special cases. Spinless fermions corre- 
spond to N = 1: matrices have only one component and 
I = 2, T = 0. Thus there is only one interaction func- 
tion. For instance, in the spinless Hubbard model with 
nearest- neighbor interaction constant Uq", the function 
U A is expressed in terms of U$ n and some combina- 
tion of trigonometric functions, depending on the spatial 
dimensionu. 

In the electron Hubbard model with on-site interaction 
constant U§ s , the functions are U A = and U s oc Uq s . 
Switching on a constant interaction E/q" 1 between nearest 
neighbors, we come up with two independent functions 
U and U in the Hamiltonian. 

The expressions (|^) and (^) for the action are ade- 
quate for a microscopic, 'exact' formulation of the prob- 
lem. The functions U A and U s may incorporate the 
microscopic interaction of our choice: Coulomb, on-site 
repulsion, and so on. The integration in Eq. (^|) is then 
carried over all available phase space (the Brillouin zone) , 
with the constraint of momentum conservation. In prin- 
ciple, working with the microscopic Hamiltonian allows 
a description of physical processes at all energy scales, 
up to atomic energies. However, such an exact solution 
is impossible in practice. We will consider the action 
defined in Eqs.(^,Q) as a low-energy effective actioifflau, 
describing correctly only those physical processes occur- 
ring at an energy scale Ao much smaller than the Fermi 
energy: 



In principle, such a low-energy effective action could be 
obtained from the microscopic action by integrating out 
(in the functional sense) the degrees of freedom associ- 
ated with the momenta lying out of a band of width 2Ao 
around the Fermi surface. 

In general, the low-energy effective action So + is 
written as an expansion in a set of Grassmann fields, with 
the requirements that the symmetries of the microscopic 
action be satisfied (for details, see Refs.[7],|). We assume 
that the effective action has the form given in Eqs (||) and 
(^) , except that the integration in Eq. (||) is restricted to 
the vicinity of the Fermi surface. The parameters of this 
action, like one-particle energy e(K) and the coupling 
functions U and U , do not coincide with those of the 
microscopic action in general. We ignore the possibility 
of interaction terms involving derivatives, or more powers 
of ip, because such terms are irrelevant at tree- level (see 
below) . 

In this work we will restrict ourselves to the study 
of low-density electron systems, so we assume that the 
Fermi surfaces have rotational symmetry (i.e. circular 
or spherical). The low- lying one-particle excitations can 
then be linearized near the Fermi energy in a simple fash- 



e(K) - p w v F {K - K F ) = v F k 



(14) 



wherein the terms of order k 2 in this expansion are irrel- 
evant. The momentum k and the temperature T are re- 
stricted by the inequalities \k\ < Ao and T <C v F k$ along 
with the condition (p^|). The bare one-particle Green's 
function for the free part of action Sq is 

(V>a(l)^/j(2)) - Go(K 1 , Wl )^ d+1 )(l - 2)800 (15) 
Gq 1 (Ki, lui) = iuji + fi — e(Ki) ~ iu)\ — v F k\ (16) 

In the (i-dimensional integration measure only the rele- 
vant term is kept: 



f ° I (K F + k) d ~ l dkd£l d 



K d-i 

F 



-An J SI. 



dkdfl. 



d 



(17) 



The Matsubara frequencies are allowed to run over all 
available values. We presume that the density of particles 
in the system is kept fixed. 

The low-energy effective action (^Ji]) lends itself to a 
simple tree- level scaling analysis (dimensional analysis). 
It is invariant under the rescaling u' n = sui n , k' = sk 
and ijj'{k' ,u>'„) — s~ 3 ^ 2 tfj(k,uj n ). Higher orders in the ex- 
pansion (|l4j) or in the measure (17) are then irrelevant 
(in the RG sense) at tree-level, and so are more com- 
plicated interaction monomials. This justifies the simple 
form (§J§). 
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According to the RG analysis of Shankara, performed 
at zero temperature, this effective action describes the 
Fermi liquid phase of repelling spinless fermions and 
the BCS superconducting phase of attracting spinless 
fermions. We will study here the RG fixed points for 
interacting A-flavor fermions expanding this technique 
for the finite temperature case. 



III. THE RG EQUATIONS IN TWO DIMENSIONS 

A. The coupling functions 

In this section we will restrict the analysis to the two- 
dimensional case, with a circular Fermi surface. As fol- 
lows from tree-level scaling, the frequency dependence 
of the coupling function U is irrelevant and is discarded 
from now on. 

Each vector Ki in the effective action lies in a narrow 
shell of thickness 2Ao around the Fermi surface. We may 
write it as K; = K F +ki where the vector k^ is orthogonal 
to the Fermi surface and small, in the sense that |k;| < 
A <C Kp. The coupling function is then written as 

jjA,S = ^A,S (K l i) K 2 . K 3^ K 4 | ki) k2 . k3) k4) (18) 

If we restrict ourselves to the case of non-singular 
interactions ,Ej the k-dependence of the coupling func- 
tion turns out to be irrelevant. Moreover, the depen- 
dence of the coupling function on the Fermi surface mo- 
menta is constrained by crystal momentum conserva- 
tion: K F + Kp = Kp + Kp (we assumed that Umklapp 
processes are inoperative, which is justified for a circular 
Fermi surface at low-filling) . 

A geometrical analysis (see Ref. ||) leads to the follow- 
ing conclusion: In the case of a circular Fermi surface, 



the constraint that the momenta K F lie on the Fermi 
surface and the momentum conservation law allow the 
existence of two pairs of independent coupling functions, 
as described below. 

Case 1. If Kp — K F , we define a dimensionless 
coupling function $ of the momenta and spins: 



$(K F , K F ) = -v F U${K F , K|; K F , K%) (19) 



where v F = 2S d Kp- 1 /(2 7r) vf is the density of electron 
states at the Fermi level and Sd is the area of the d- 
dimensional unit sphere. The hat (*) means that $ is 
an operator in flavor space. Each vector K F may be 
specified by a plane polar angle 9 a ; because of rotation 
invariance, the function $ may only depend on the rela- 
tive angle 9\ — 9 2 between Kp and Kp: 



-0 2 ) - * A (0i-e 2 )tf 



S ' a i-6 2 )T^ (20) 



It should be pointed out that the function <& A (9) is not 
an antisymmetric function of its argument. It follows 



from the symmetry properties of the coupling function 
(Eqs. (§,|l|l|)) that 

^ s (9 1 -9 2 ) = <S> A ' S {6 2 -Qx) (21) 

The only remnant of the antisymmetry of U A is the con- 
dition $ A (0) = 0. It is easy to check that an exchange 
of incoming or outgoing momenta does not require any 
new independent function. 

Case 2. The two outgoing momenta are opposite: 
Kp (likewise for the two incoming momenta). 
We introduce another pair of dimensionless functions: 

1 



K 2 



V{Kp,K F ) 



:VfU$(K\ 



KIiKL-K? 



p,— xv F ,xv F , xv F ^ (22a) 

V S {9^)T^ (22b) 

wherein #13 = 0\ — #3. From the symmetry properties of 
the coupling function U under exchange of momenta, we 
find 



V a (9±tt) 



■±ir) 



-V A (d) 
V s (6) 



(23a) 
(23b) 



There is an overlap in the definitions of $ and V, since 
3>(±7r) = V'(O) or, more explicitly 

$ A (±vr) = V A (0) = -V a (±tt) 



$ s (±tt) = V s (0) = V s (±tt) 



(24a) 
(24b) 

Because of the 27r-periodicity of the functions (p" 4 * 5 ', V A ' S 
and of the symmetry properties ( pl|) and (|2^) , these func- 
tions take all their possible values in the following 'min- 
imal' domains: 9 e [0, tt) for & A ' S (9) and 9 6 [0, tt] for 
V A ' S (9). The angle 9 = 7T is excluded from the domain 
of & A ' S in order to avoid any ambiguity following from 

©■ 

In the rotationally invariant case it is convenient to 
expand the coupling functions in Fourier series: 



X(0) 



E 

l= — oc 



X,= f' y X(9) 
Jo 2vr 



no 



(25a) 



(25b) 



where X stands for the set of all coupling functions: X = 
$ s , V A , V s }. In terms of Fourier components, the 

a A S A S 

symmetry properties of <& become $ ; ' = $ _j together 
with 



E 



$f = 



(26) 



Those of V become V A — (I even) and Vf = (I 
odd). The symmetry properties of the coupling functions 
expressed above are exact, for they are consequences of 
the Pauli principle. 

Let us finally point our that the functions Q A and V A 
coincide with those introduced by ShankarEI in the spin- 
less case, while $ s and V s are brought about by the 
introduction of spin. 
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B. RG equations in the Landau channel 

Consider the two-particle Green's function: 

G 2 (l,2;3,4) = -(^1)^/3(2)^(3)^(4)) (27) 

The vertex function 

r(l, 2; 3, 4) associated to G2(l, 2; 3, 4) is constructed by 
considering only connected one-particle-irreducible (1PI) 
diagrams with amputated external legs. In the first-order 
approximation, fW (1, 2; 3, 4) = (7(1, 2; 3, 4). 
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-111 

2 , 






(ZS) 



/' 




(ZS') 



(BCS) 




(ZS2) 




irrelevant 
diagrams 



(BCS2) 



FIG. 1. Perturbative diagrammatic expansion for the total 
vertex in terms of the bare interaction. At the two-loop level, 
only the diagrams employed in the Sec. HID are depicted. 

Let us define a flow parameter t such that A(i) = 
Aoe _t . We will follow the flow of the vertex function from 
t = to t = oo. We introduce the dimensionless tem- 
perature /9fl(t) = ^vp0A(t). Three Feynman diagrams 
contribute to the Gell-Mann-Low (GML) function of the 
vertex function at the one-loop level (see Fig. 1): zero 
sound (ZS), Peierls (ZS'), and Cooper pairing (BCS). 
The coupling functions <f> and V can be extracted from 
the vertex T by proper choice of the external momenta 
lying on the Fermi surface and by putting the external 
frequencies equal to the minimal Matsubara frequency 
cumin = ttT '. Summations over flavor indices and Mat- 
subara frequencies can be_.done easily using Eq. jic| ) 
and standard techniquesa'til for working with tempera- 
ture Green's functions. The relevant contributions to 
the GML functions of $ and V come only from those 
diagrams in which all the momentum phase space is avail- 
able for integration. For a detailed discussion of this last 
issue, see Ref. |[ „ 

To our knowledge, Shankar was the first to argueu 
that the RG fixed point in the Landau channel yields 



the Landau /-function of the Fermi liquid. At zero tem- 
perature, Shankar obtained an RG-invariant function $> A 
and identified it with the Landau function. Our finite- 
temperature analysis reveals a more subtle situation. 

As is known from the Landau FL theorytl, there are 
two limits of the four-point vertex function with all four 
momenta lying on the Fermi surface, depending on the 
way the limit of zero momentum- and energy-transfer 
is taken. One of these limits, which does not possess 
the symmetry (26) of the forward scattering amplitude , 
yields the Landau function. The other one, which pre- 
serves the antisymmetry of the vertex, gives the total for- 
ward scattering amplitude. For details see, Jar instance, 
the especially elucidative paper by MerminJl^l 

To clarify this point in our approach, let us consider 
the perturbative contribution of the ZS diagram to the 
vertex function f(l, 2; 3, 4) at small momentum transfer 
q = K3 — Ki and small energy transfer iQ 13 = i{u!i —UJ3). 
As discussed after Eq. (p!s|), we neglect the difference K l — 
K F in the effective interaction. For this contribution we 
have: 



i zs (K F ,K F , q ,n 13 ) = J ^f/;"(4,K f ;4,K f ) 

xU v J(K F , K|; K F , K F )L ZS (K, q, n 13 ) (28) 
where 



L zs (K,d,n 13 ) 

x cosh - ' 



1 smh(/3v F K ■ q/2K F ) 

2 £Oi3 — «j?K ■ q/K F 



-/3v F k 



cosh 



I3vf 
2K F 



K F • (k + q) 



(29) 



is the result of the summation of the product of two 
Green's functions over Matsubara frequencies inside the 
loop. As follows from Eq. (p9|), Lzs has two limits: In 
the 'unphysical' limit (q — > 0, followed by SI13 = 0) Lzs 
vanishes. In the 'physical' limit = 0, followed by 

q -» 0), L zs = -i/3cosh- 2 (^ F fc). 

All the phase space is available for integration in the 
r.h.s. of Eq. (|28|), so the ZS diagram gives a relevant 
contribution to the total vertex at any angle 9 between 
JCp and K F . A similar analysis for the ZS' diagram 
shows, that because of the phase space restrictions, it is 
relevant only when 9 — ► 0, and its contribution does not 
depend on the way how zero energy- and momentum- 
transfer limit is taken: 

lim X ZS '(0,0) = J -^U^(K F ,K F ;K F ,K F ) 

xU^(K F ,K F ;K F ,K F )^pcosh~ 2 (^l3v F k) (30) 

In the following we will call &(9,t) the renormalized di- 
mensionless vertex in the 'unphysical' limit, i.e. 



mt) = h F 



lim r(l,2;3,4) 

q^0 



(31) 



Oi 3 =0 
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This notation is justified since the vertex function coin- 
cides with the effective interaction $ in that limit. On 
the other hand, we define the renormalized dimensionless 
vertex in the 'physical' limit as 



f(l,2;3,4) 









(32) 



We will use the components of the vertex T A ' S defined 
in the same fashion as in Eq. (|TT|). 

Note that among one-loop contributions to the GML 
function of the vertex F with arbitrary incoming and 
outgoing external momenta and frequencies, the BCS 
channel possesses the total symmetry of the interac- 
tion (Eqs. (|l2|)), while the ZS (or ZS') channel sepa- 
rately do not: only their combined contribution is an- 
tisymmetric under exchange of incoming (or outgoing) 
particles. Here, the diagram ZS' contributes to the 
GML function of &(9,t) when 9^0, and it generates 
the flow for <fr A (9,t) at 9 = 0, breaking the condition 
<I>^(0) — 0, while this ZS' contribution for the physi- 
cal vertex cancels dT A (9,t)/dt at 9 = 0. This makes 
the derivatives d$ A (9,t)/dt and dT A (9,t)/dt discontin- 
uous at 9 = 0. In reality, this feature is smeared over 
some angle 9 sm ~ A /Kp, which is small (9 sm <C 1) be- 
cause of the condition (13). In order to obtain a smooth 
GML functions for ${9,t) and t(9,t), we adopt the fol- 
lowing procedure: First, we define <t(9,t) and T(6,t) in 
the patch 9 £ [9o,ir), where 9 sm < 9q <C 1. Then we 
derive the RG equations in this domain, and finally we 
analytically continue <£(#, t) and T(9, t) up to 9o = 0. 

The calculations are straightforward, and we end up 
with the following RG equations: 



9$ 



A,S 



dt 



= o 



or 



A,S 



iA,S 



dt 



cosh 2 /3r 



dt 



Here we introduced the following notation: 

N — 1 



P A = 



N+l 



[r 



Ai2 
I \ 



(N- 
■(N 



A-rS 



i)rfr 



Ft 



Si 2 



1 



irfrf 



N - 3 



[r 



Si 2 
/ \ 



(33a) 
(33b) 
(33c) 

(34a) 
(34b) 



As follows from Eqs. (|33|a), the functions $ A,S (9, t) are 
RG invariant. The effect of the singular contribution 
to the RG flow from the ZS' diagram is not included in 
Eq. (p3|a) , but is equivalent to relaxing the symmetry con- 
dition Q A (0) = for the bare coupling function & A (9). 
It is convenient to express the RG invariants $> A ' S in 
terms of 'charge' (F) and 'flavor' (G) functions, defined 
as follows: 



F= (N - 1)$ 5 - (N + l)& 
G = - $ A 



(35a) 
(35b) 



The physical meaning of these functions becomes clear in 
the case of real electron spin (N = 2), when theycoincide 
with the components of the Landau /-functionEJ: 



/i2,3 4 (0) - F^(9)S 13 8 24 + G (L) {9) 



0"13 • CT24 



(36) 



In order to decouple the Eqs. (p3|b) we introduce new 
components of the renormalized vertex T: 



A = (N - 1)F 5 - (N + l)T J 

b = -r s - r A 



(37a) 
(37b) 



Their RG fixed points in the case N = 2 coincide with 
the components of the tatal scattering amplitude in the 
notations of FL theory.E3 (Notice, that the components 
r s and — T A are respectively called singlet and triplet 
amplitudes in Ref. |l6|). Notice also the symmetries Ai 
A-i and Bi = B-i. To simplify the system of Eqs. ((33|t 
and (p3c) we use an auxiliary RG parameter 



r = tanh/?R , r £ [0,Tq] 



(38) 



wherein tq = tanh/?o < 1, and /3q = ^vf(3A.q. 

Expressed in terms of these new variables, the RG 
equations become much simpler: 



dr ~ Al 
9Bl B 2 



(39a) 
(39b) 



From Eqs. (|39|) we obtain the following stable fixed points 
(i.e. the solutions at t = oo) : 



A*i = 



F, 



1 + FiT 
Gi 



1 ~ 1 + G ; r 



(40a) 
(40b) 



where we expressed the bare values of the vertex via pa- 
rameters of the effective interaction (the FL parameters) , 
Gi. (In this work asterisks and 



.4 



(0) 



F t , B 



(0) 



zeros respectively denote fixed points and bare values.) 
From Eq. (Q) it follows that, if the parameters of inter- 
action satisfy the conditions 



{f,,G/} > -1 



V I 



(41) 



the system will remain stable in the Landau channel 
at any temperature. The fixed point (^o|) at T = 
(to = 1) gives the same solution for the forward scat- 
tering vertex as the BethpJ3alpeter equation in the zero- 
temperature technique 3V3 Notice that, because of the 
condition T <^ vfAq imposed on the low-energy effective 
action, we can set To = 1 for all practical purposes. Any 
attempt to extract a critical temperature from Eq. |40j 
by looking at the violation of condition ([H]) for the ab- 
sence of poles in the vertex function, would involve ex- 
ponentially small corrections to zero temperature, thus 
exceeding the accuracy of the present approach. 



G 



To take into account the RG flow from the ZS' diagram, 
which preserves the symmetry condition T A (8 = 0,t = 
oo) = 0, we impose the constraint 



[(N-l)At + Bt) = 



(42) 



on the fixed point of the vertex. This specific form ( M^ ) 
of the condition for the forward scattering amplitu<jkjs-|is 
usually called the amplitude sum rule in FL theory.EH13 



C. The RG equations in the BCS channel 

The only relevant contribution to the GML function 
for V(8,t) comes from the BCS diagram. The fact that 
the BCS diagram preserves by itself the symmetry (|23| ) 
of the function V makes life simpler. The RG equations 
are derived easily: 



dt 



tanh/?ij (V t A ) 2 



— _ -tanhfo (Vf) 2 



(I odd) (43a) 
(/ even) (43b) 
) coincides with Shankar's 



At zetp temperature, Eq. (ktj 
resulta for spinless fermionsl 

The RG equations ( f43| ) may be solved easily, and their 
fixed points are: 



V, 



A(*) _ 



V, 



A{0) 



r °°tanh(/3 e-*)* 



V, 



s(*) _ 



V, 



S(0) 



l + V, s<p) / o °°tanh09be-*)<ft 



(44a) 
(44b) 



The system remains stable in the V channel at any tem- 
perature if, for all harmonics, the following conditions 
are fulfilled: 



V, 

v, 



Mo) 

S(0) 



< 
> 



(45a) 
(45b) 



At low temperature (i.e., when the low-energy action ap- 
proach is supposed to work well), (3q 3> 1 and the inte- 
grals of Eqs. (|44|) can be evaluated exactly, like in the 
theory of superconductivity (cf. section 33.3 of Ref. ||): 



A(*) 



V, 



V, 



A(0) 



1 - \n{2^v F A a /irT) 



V, 



S(0) 



1 + vf (0) ln(27i; F Ao/7rr) 



(46a) 



(46b) 



wherein In 7 = 7 ks 0.577 is Euler's constant. If the 
conditions (^|) are violated, the ^-interaction becomes 
marginally relevant, and a pole appears at temperature 



Tgc — — «FAoexp(- 



2 7 



(47) 



Here S/ = max^V^^ , |V^ S( ' ' ) |}, in which only the har- 
monics violating the condition ( p5| ) are included. If 
phonons provide the underlying mechanism of attraction, 
we may identify the characteristic energy scale vf-^-o of 
the low-energy action with the Debye energy uor>- In 
the present context, the BCS theory of superconduc- 
tivity corresponds to the special case of a contact at- 
tractive interaction, for which only the zeroth harmonic 
oc Uq < does not vanish. 



D. Beyond the one-loop approximation 

The main goal of this subsection is to argue that there 
are no higher-loop corrections neither into one-particle 
Green's function, nor into the vertex. In other words, 
the one-loop RG results are exact_for this model, based 
on the low-energy effective action.El 

The phase space analysis of the diagrams contribut- 
ing to the self-energy E shows, that the only relevant 
contribution comes from the tadpole diagram at the one- 
loop level, and diagrams obtained from it by tadpole self- 
energy insertions at higher loops. This self-energy con- 
tribution is temperature-independent, and it results in a 
mere shift of chemical potential of the interacting system 
flint- As far as we retain the constant density of parti- 
cles in the system, the self-energy contribution exactly 
compensates the change of the chemical potential due to 
interaction at T — 0, i.e. //j n t(T = OWE = fi = Kp/m* , 
according to the Luttinger theoremO So, the renormal- 
ization of the one-particle Green's function due to the 
self-energy holds its linearized form ( |l6| ) near the Fermi 
surface unchanged. It supposes also the temperatures 
under consideration being small enough to neglect the 
temperature corrections to the chemical potential of the 
interacting system. 

Thus, integration over the rest of the degrees of free- 
dom inside the narrow shell around the Fermi surface 
neither generates an additional wavefunction rcnormal- 
ization (we have chosen the initial scale of the Grass- 
mann fields in the bare effective action such that Z = 1), 
nor an additional renormalization of the effective mass 
of quasiparticles (i.e. the phenomenological parameter of 
the effective action vp = Kp/m* does not change). 

To show that the RG results obtained above are ex- 
act for this model, we will demonstrate the cancellation 
of the next terms in the GML functions at the two-loop 
level. The generalization of this proof to higher loops is 
straightforward. In this subsection we use the standard 
field-theory renormalization techniquetalij with momen- 
tum integrations taken over the interval k G [—A, A], in- 
stead of the Wilson-Kadanoff (WK) scheme. Since we 
work at finite temperature, there is no need to introduce 
an infrared cutoff, because all integrals are regular near 
the point k = 0. The pictorial form of the equation for 
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the renormalized vertex in terms of the lPI-diagrams is 
presented in Fig. 1. The light squares stand for the bare 
vertices, while the dark square denotes the renormalized 
vertex. The equation for the renormalized vertex may be 
written in the symbolic short-hand form as: 

ur = u - u 2 Li + WqL 2 (48) 

From Eq. (g§) 

the GML function @gml = (dun/d\nA) UQ is found to 
be 

Pgml{ur) = -u 2 R L[ + u 3 R (L 2 - L\)' (49) 

where the primes in the r.h.s. of Eq. ( f49| ) stand for the 
partial derivative w.r.t. In A. We can use symbolic form 
( [49| ) for the calculation of the GML functions for renor- 
malized functions $(0, A), f (0, A), and V(9, A). 

The phase space analysis of the two-loop diagrams 
shows that the only relevant contribution to L2 in the 
Landau channel comes from the diagram ZS2. In the 
unphysical limit this contribution is zero, while in the 
physical limit it is exactly canceled by the term — L\ 
which in this case is the contribution from the one-loop 
ZS diagram. Therefore, Eqs. (p3|a, b) are true at two-loop 
level as well. In the BCS channel, the only relevant con- 
tribution to L 2 in the equation for the GML function of 
V(0,A) comes from the diagram BCS2, and in the same 
way it is canceled by the contribution L\ from the BCS 
diagram, leaving Eq. (|Z^) unchanged. 

It is worth noticing that the above proof of the cancel- 
lation of higher-order terms of the GML functions do not 
work in one dimension, because the phase space restric- 
tions in ID are not strong enough to remove the contri- 
butions from all higher-loop diagrams. At the two- loop 
level, the term Li in Eq. ( |49| ) includes the contributions 
from all 1PI vertex diagrams along with the self-energy 
'sunrise' diagram contribution, responsible for the wave- 
function renormalization. Thus, in ID there is no can- 
cellation in the GML functions due to renormalization. 



IV. THE RG EQUATIONS IN THREE 
DIMENSIONS 

In order to obtain the marginal coupling functions for 
a spherical Fermi surface in 3D, we follow the same pro- 
cedure as in Sec. III. There are two marginal coupling 
functions: 

1. The function l>. It is defined by Eq. (p"9"|). It couples 
a pair of fermions with momenta (K-,K|7) to another 
pair with momenta (K^,K^), lying on the cone swept 
by a rotation of (K-, Kf.) around the vector + K.p. 
Let 0i2 be the angle between two vectors and Kl, 
and 0i3 the angle between the two planes defined by the 
pairs (K F ,K 2 F ) and (K^K^). The function $ may be 
written as 

$ = $>i 2 , <h s )I°i + $ s (0 12 , 0i 3 )Tf (50) 



wherein the functions $> A < S have the following symmetry 
properties: 

$ a (012,0i 3 ±tt) = -$ A (0i 2 ,0i 3 ) (51a) 
$ s (0i 2 ,0 13 ±7r) = $ s (0i2,</»i3) (51b) 

The minimal domain of definition allowing to recover all 
values of these functions is 0, £ [0, vr). 

2. The function V . It is defined by Eq. (p2|a), with 
the symmetry properties (p3[). 

The phase space analysis shows that the RG flow is 
generated only in the forward scattering amplitude chan- 
nel, i.e., in the 'physical' limit of the vertex T(0). All 
other coupling functions are RG invariant and never mix 
up with the = channel. Like in the 2D case, the rel- 
evant contributions to the GML functions of T(9,t) and 
V(6, t) come respectively from the diagrams ZS and BCS. 

The spherical symmetry of the Fermi surface makes 
convenient an expansion of these functions into Legendre 
polynomials. Using again the short-hand notation X for 
the set of coupling functions (cf Eq. (p5|)), this expansion 
is 

00 

X(0) = ( 2l + !) X ' A(cos0) (52a) 
1=0 

X; = - / d(cos0) X(0) F ; (cos0) (52b) 
2 J-i 

Here, in the set X(0), <S> A < S (6) stands for $ A ' S (0, = 0). 
All the calculations are conducted in the same way as in 
the previous section. 

The symmetry properties and the RG equations of V 
are exactly the same as in two dimensions. The only 
technical modification in the Landau channel is the sum 
rule for partial amplitudes (cf. Eq. (E^) 

00 

£>l + l)pr-l)A? + B,*]=0 (53) 
;=o 

The rest of the formulas and all the conclusions made in 
the 2D case may be repeated word for word in the 3D 
case. 



V. RESPONSE FUNCTIONS 

In this section we apply our RG technique to the study 
of response functions. (For the definitions of these func- 
tions and their relationship with physical observables, 
see, for instance, Ref. |ll|.) All the results presented be- 
low are valid for spatial dimensions d = 2,3. In order 
to calculate the compressibility, we consider the density 
response function k(Q), defined as follows: 

<^ +1 >(0MQ) = (/5(Q)p t (Q)} (54) 
wherein the density operator p(Q) is 
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P(Q) = [ ^a(l)^a(l + Q) 
J(l) 



(55) 



and p(Q) = p(Q) — (p(Q)) stands for the density fluc- 
tuation. Here Q = (q, f2„), and il„ is the Matsubara 
frequency. The zeroth component of this function k(Q) 
taken in the physical limit gives us the derivative of the 
particle concentration n w.r.t. the chemical potential, 
i.e. 



k = lim 



«(Q) 






o„=o 



dn 



(56) 



The RG flow is simplest when obtained in the field- 
theory renormalization scheme, like in the Sec. HID. 
Notice, that in this scheme the RG parameter r (cf. 
Eq. $3^)) runs from t = to r = t « 1 (the fixed 
point), and the r.h.s. of Eqs. ( |39| ) changes its sign. The 
first two terms of the perturbative expansion for n give 
us (cf. Eq. (137" 



(57) 



wherein kq = ^vpT is the contribution from the free 
part of the effective action. According to the results of 
the Sec. HID, to obtain the equation for the RG flow 
of k(t), exact in this model, it is sufficient to take into 
account one-loop renormalization of the physical vertex. 
Introducing the auxiliary function k(t) = k/kq, and us- 
ing Eqs. (|39|), we get the RG equation: 



<91nft(r) 
dr 



= -Mr) 



which yields dn/dfi as its fixed point: 

, = N V F 

2 1 + F 



(58) 



(59) 



From the thermodynamic formula for the compressibil- 
ity K = -(l/V)(dV/dP) = (l/n 2 )(dn/dfi) (see, for in- 
stance, Ref. M), we easily recover the result for the elec- 
tron Fermi liquid: 



K = 



1 Vp 



(60) 



To find the spin susceptibility in our SU(N) formalism, 
we consider the flavor response function x(Q)> defined as 
follows 



S (d+1) (0)x(Q) 



f 



N 2 



T (S„(Q)5t(Q)) 



wherein the flavor density operator 



S a (Q) = 



(1) 



A»^ Q (1)^(1 + Q) 



(61) 



(62) 



and g is the gyromagnetic ratio. The uniform suscep- 
tibility x i s given by the physical limit x(Q — > 0) (cf. 



Eq. d56|)). Defining x(Q) m the fashion (^lj), we used 
the fact that in the paramagnetic state the response is 
the same along all the (./V 2 — 1) directions a. The rest of 
the calculations is carried out in the same way as above. 
For the auxiliary function x(t) = x/xo (xo = \g lv FT is 
the contribution in the absence of interaction) we obtain 
the RG equation: 



dhix(T) 
dr 



-B (t) 



Again, we recover the FL theory result 

Vp 



E 



X. 



4 y 1 



Go 



(63) 



(64) 



Notice, that the stability conditions for the fixed points 
(p9|) and (|64|) a re just a special case (I = 0) of the general 



conditions ([41 






+ 



FIG. 2. Diagrammatic form of the recursion relations for 
the susceptibility x(t) (circle) and the vertex z(t) (triangle). 
The hatched symbols stand for the functions at t+dt, whereas 
light ones stand for functions at t. The black square stands 
for the vertex T(t). Integration over loop momenta is carried 
out in an infinitesimal shell dA(t). 

It is interesting to show how the results (|9|) and Q) 
can be obtained in the Wilson-Kadanoff scheme.LJ Below 
we specialize to the electron (N — 2) spin susceptibility, 
but the compressibility may be computed in the same 
fashion. We add to the effective action a source field 
/i(Q), conjugated to the z-component of the spin den- 
sity (p2| ). For details on this approach, see, for instance, 
Ref. ^The successive integration over momentum modes 
generates a vertex correction z[t) to the source term 
along with higher order terms in the external field h. 
The effective action, at some intermediate value of t (the 
flow parameter), takes the form 

S[i>, h, t] = S[il>, 0,t}+ f z(t, Q) [h(Q)S 3 (Q, t) + h.c] 

J(Q) 

+ [ X (t,Q)h*(Q)h(-Q) + 0(h 3 ) (65) 

J(Q) 
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In linear response theory it is sufficient to keep track of 
terms up to second order in h. The recursion relations for 
the vertex z(t) and the susceptibility x(t) are illustrated 
on Fig. 2. In the physical limit Q — * 0, we obtain the 
following pair of RG equations: 



d X 1 2 2 / \ 

- = -~g v F z (r) 



g = z(r)S (r) (66) 



with the initial conditions z(tq) = 1 and x( T o) = 0- Us- 
ing the solution for Bq(t) which follows from (j39|b), we 
can easily solve this system. The uniform susceptibility, 
which is the fixed point of these equations, is again given 
by Eq. (|64|). Notice that one cannot simply compute the 
susceptibility (or any other response function) from the 
usual diagrammatic expression at the fixed point action: 
one would obtain zero at any non-zero temperature. 

In order to obtain the heat capacity in this model, one 
does not need to perform any complicated calculation. 
As explained in Sec. HID, tracing over fermion modes in 
the vicinity of the Fermi surface in d — 2 and 3 does not 
modify the one-particle Green's function ([i"6|). Therefore, 
the only relevant part of the free energy-, expressed as a 
functional of the total Green's functiontll is the same as 
for the effective action (^) without the interacting part. 
It simply gives the Fermi liquid result C = ^ir 2 NvpT. 



VI. CONCLUSION 

The RG approach to interacting fermions has been ex- 
tended to the case of non-zero temperature and spin. We 
studied a model with 5?7(A)-mvariant short-range effec- 
tive interaction and rotationally invariant Fermi surface. 
In general, the RG approach turns out to be equivalent 
to Landau's mean-field treatment of the Fermi liquid. 



The first set of conditions (45) for the bare couplings of 
the low-energy action is nothing but Landau's theorem 
for the stability of Fermi liquids .a&ainst Cooper pairing 
at arbitrary angular momentumJIHI Likewise, conditions 
( pd"| ) for the system to be stable up to zero temperature 
are Pomeranchuk's conditions for the components of the 
Landau interaction function. Only if both conditions ( |45| ) 
and ([y]) are fulfilled does the system reach zero tem- 
perature in the Fermi liquid regime. One advantage of 
the finite-temperature formalism is to reveal the possible 
breakdown of the FL picture below some critical temper- 
ature T Cl whereas it remains applicable above T c . 

The distinction between physical and unphysical limits 
as momentum and frequency go to zero appears natu- 
rally in this finite temperature formalism, and differenti- 
ates the Landau function from the scattering vertex. The 
former is RG invariant, whereas the latter flows towards 
a fixed point (^0|), which reproduces the well-known re- 
lationship between components of the Landau function 
and the forward scattering amplitudes. Working in the 
finite-temperature formalism made explicit the equiva- 
lence of the results obtained in the Wilson-Kadanoff and 



field-theory schemes of renormalization. Notice that, in 
the Wilson-Kadanoff scheme at zero temperature, the 
RG flow of the-.physical scattering vertex might be eas- 
ily overlooked^ since special care must be taken when 
the shrinking cut-off crosses the scale given by non-zero 
transferred momentum. 

In general, temperature plays the role of a 'natural 
regulator' in this theory. Firstly, it eliminates the need 
for an infrared cut-off in the BCS channel. Secondly, and 
more important, it makes the Gell-Mann-Low functions 
smooth in the Landau channel.l Thus, the fixed points 
are smooth functions in the T — > limit. 

Response functions, like the spin susceptibility or the 
compressibility, are naturally defined in the physical 
limit, and thus their RG flow involves the physical ver- 
tex r(t), which then plays the role of a 'running' effective 
interaction. This makes superfluous the use of any ad- 
ditional perturbative calculation!] once the fixed point is 
reached. 
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* The reader might be slightly confused by the sign of 
the r.h.s. of Eq. (ft^a) when comparing our result with 
that of Shankar for the spinless Hubbard model. In fact, 
this sign is just a consequence of the way we define the 
coupling function U in this work. Indeed, let us con- 
sider the Hubbard Hamiltonian on a d-dimensional hy- 
percubic lattice (with lattice constant a) at low filling, 
with on-site (Uq 3 ) and nearest-neighbor (Uq 71 ) interac- 
tions. Both constants are positive for a repulsive inter- 
action, and negative for attraction. Fourier-transforming 
and symmetrizing the interaction, we end up with the fol- 
lowing coupling functions of the microscopic Hamiltonian: 
(7o A (Ki,K 2 ;K3,K4) ~ -\a 2 US n ■ (Ki - K 2 )(K 3 - K 4 ) 
and (7 s (Ki, K 2 ; K 3 , K 4 ) ~ U% 3 +U$ n . Here, only the low- 
est order terms in the expansion of trigonometric functions 
were kept. Then we can use this result as the lowest-order 
approximation for the bare coupling functions of Sint in 
the definitions @ and (||). 
** For the following discussion we are indebted in great part 
to N. Dupuis. 

1 Notice that, in Eq. (||b), lim T _ (l/4T) cosh~ 2 (:r/2T) = 
8(x). 
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